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We discuss the production of 26 Al and 60 Fe in two massive star models, namely a 25 
and a 80 M , of initial solar metallicity. 

O 

1. Introduction 

o , 

The radioactive 26 Al and 60 Fe isotopes decay into their daughter stable nuclei with 
lifetimes of the order of a few ~ 10 5 yr. These times are short enough, compared to the 
typical timescales of Galaxy evolution, that the observations of these isotopes clearly iden- 
tify the locus of the present nuclesynthetic enrichment through the Galaxy. Observations 
of the 7 — ray emission from the radioactive 26 Al at 1.809 MeV have been reported by 
several missions over the last 15 years and show that 26 Al is mostly confined in the disk of 
the Galaxy with some prominent areas of emission corresponding to the Vela and Cygnus 
regions (Diehl & Timmes 1998). Many sources can be responsible for the production of 
26 Al in the Galaxy, among them Core Collapse Supernovae, WR stars, Novae and AGB 
stars. Unfortunately the observations do not provide a strong constraint on the dominant 
source of production yet. There is another measurable quantity, however, that could help 
to distinguish between the various competing sources. In particular 7 — rays from the 
radioactive 60 Fe may be a good discriminant for the origin of Galactic 26 Al. Indeed, core 
collapse supernovae should produce comparable amounts of both 26 Al and 60 Fe while all 
the other candidates produce a much smaller amount of 60 Fe compared to 26 Al. 

A few years ago we started a series of papers devoted to the study of the presupernova 
evolution and explosive nucleosynthesis of massive stars, in the range 13-35 M , as a 
function of initial mass and metallicity (Chieffi et al. 1998, Limongi et al. 2000, Chieffi 
& Limongi 2004). No mass loss was included into those computations. At present we are 
addressing the presupernova evolution and explosion of massive stars with mass loss in 
a more extended range of mass, namely between 13 and 120 M (Chieffi & Limongi in 
preparation). In this paper we preliminarily discuss the production of 26 Al and 60 Fe in 
two selected massive star models, namely a 25 and a 80 M , of initial solar metallicity. 
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Figure 1. Upper Panel: internal profile of 26 Al and 60 Fe (thick solid lines) for the 25 M Q 
model before the explosion. The dotted lines refer to the most abundant isotopes. The 
y-axis of these lines is not shown and ranges between and 1. Lower Panel: internal 
profile of 26 A1, 60 Fe together with their cumulative values (thick solid lines) for the 25 M 
model after the explosion. Also in this case the most abundant isotopes are shown. The 
y-axis of 26 Al and 60 Fe is not shown and ranges between and 5 ■ 10 -4 . 



2. Production of 26 Al and 60 Fe in massive stars 

In stars, 26 Al is produced by proton capture on 25 Mg and destroyed by (3 + decay into 
26 Mg and by (n,p), (n,a) and (p, 7) reactions (Clayton & Leising 1987). In massive stars 
(M > 11 M ) 26 A1 can be produced by (1) hydrostatic H burning, (2) hydrostatic C 
and Ne burning and (3) explosive C and Ne burning. Any 26 A1 produced by these stars 
is ejected by both stellar wind and explosion in different proportions depending on the 
initial mass. 

60 Fe is produced by the sequence 58 Fe(n, 7) 59 Fe(n, 7) 60 Fe. This sequence is mainly pow- 
ered by the 22 Ne(a, n) 25 Mg reaction, that is efficient during (1) hydrostatic He burning, 
(2) hydrostatic C burning and (3) explosive Ne burning. 

In this paper we discuss the production of 26 Al and 60 Fe in two massive star models, 
namely a 25 and a 80 M , of initial solar metallicity. These models have been computed 
with the latest version of the FRANEC code (Chieffi & Limongi 2004). Mass loss have 
been included following the prescriptions of (1) Vink et al. (2000, 2001) for the blue 
supergiant phase (T e fj > 12000 K), (2) de Jager et al. (1988) for the red supergiant 
phase (T eff < 12000 K) and (3) Nugis & Lamers (2000) during the Wolf-Rayet phase. 
The explosion has been computed adopting a PPM hydro code in the framework of the 
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piston method. The initial velocity of the piston is set in order to eject all the matter 
above the iron core. As a consequence the final kinetic energy is different depending on 
the progenitor mass. 

Figures 1 and 2 show the internal profiles of both 26 Al and 60 Fe for the 25 and a 80 M 
models respectively. The upper panels refer to the pre explosive abundances while the 
lower panels refer to the final explosive ones. 
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Figure 2. Same as Figure 1 but for the 80 M model. The y-axis of 26 Al and 60 Fe in the 
lower panel is not shown and ranges between and 6 • 10~ 4 . 



In the 25 M star the bulk of 26 A1 is produced by the Ne shell during the presupernova 
evolution; then it is completely destroyed by the explosion and reproduced, farther out 
in mass, by the explosive Ne and C burning. Hence, in this case, the final yield of 
26 Al (1.69 ■ 10 -4 M ) is dominated by the explosive nucleosynthesis. On the contrary 
the 60 Fe produced by the C convective shell during the presupernova evolution is only 
marginally affected by the explosive Ne and C burning. Hence the total final yield of 
60 Fe (1.15 ■ 10~ 4 M ) is due, by more than 70%, to its hydrostatic production. Let us 
eventually note that no 60 Fe is synthesized in a sizeable amount within the He convective 
shell. 

In the 80 M star the production of 26 Al is more complicated than in the 25 M . 
Indeed 26 A1 is initially produced during the presupernova evolution by the central H 
burning. Roughly half of this 26 A1 is lost by stellar winds during the Wolf-Rayet phase, 
while the remaining is burnt during the following He burning. Another episode of 26 Al 
production occurs during the subsequent C and Ne shell burning phases so that the final 
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profile of 26 A1, that is present within the star at the moment of the explosion, is shown 
in the upper panel of Figure 2. The effect of the explosion is to completely destroy the 
26 Al produced by the hydrostatic Ne burning shell and to reproduce it farther out in 
mass. The final total yield of 26 Al (5.13 • KT 4 M ) is due to the sum of the 26 Al present 
in the stellar wind (1.11 • 10~ 4 M ), the one that survived the explosive nucleosynthesis 
(1.29 • 10~ 4 M Q ) and the one due to the explosive burning (2.73 • 10~ 4 M Q ). At variance 
with the 25 M model, the He convective shell dominates the final nucleosynthesis of 60 Fe 
(see the upper panel of Figure 2). Such an occurrence is due to the fact that, in this case, 
the He convective shell is exposed to a temperature high enough that the 22 Ne(a, n) 25 Mg 
reaction becomes extremely efficient. Such a high temperature obviously induces also a 
large He consumption (~ 0.195 in mass fraction) within the He convective shell. The 
final total yield of 60 Fe (9.09 • 10~ 4 M ) is due to the sum of the 60 Fe produced by the He 
convective shell (7.99-10 -4 M ), the one produced by the C convective shell (6.4-10~ 5 M ) 
and the one produced by explosive Ne and C burning (4.6 • 10~ 5 M ). No 60 Fe is present 
within the stellar wind. 

3. Summary and conclusions 

We have discussed the production of 26 Al and 60 Fe in two massive star models, namely 
a 25 and a 80 M , of initial solar metallicity. We have shown that in the 25 M model, 
the final yield of 26 A1 is dominated by explosive C and Ne burning while 60 Fe mainly 
comes from the convective C shell (~ 70%) and to a lesser extent (~ 30%) by explosive 
C and Ne burning. On the contrary, in the 80 M model 26 Al comes from: stellar wind 
(~ 20%), convective C shell (~ 25%) and explosive C and Ne burning (~ 55%). The 
final yield of 60 Fe is dominated by the He convective shell by more than 85%. No 60 Fe 
is present in the stellar wind. We explored the dependence of the final 26 A1 yield on the 
explosion energy and we found that it remains essentially constant up to, at least, final 
kinetic energies of the order of a few foes. 

As a final comment let us note that the 60 Fe/ 26 Al ratio is 0.68 and 1.77 for the 25 
and the 80 M models respectively. These ratios are both well above the upper limit 
(~ 0.20) reported by the 7-ray instruments presently in Space. A detailed discussion of 
the theoretical steady state 60 Fe/ 26 Al prediction, in the light of the present uncertainties 
in the computation of the stellar models, will be given in a forthcoming paper. 
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